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It was recently shown that, for solving NP-complete problems, adiabatic paths always exist with- 
out finite-order perturbative crossings between local and global minima, which could lead to anti- 
crossings with exponentially small energy gaps if present. However, it was not shown whether such 
a path could be found easily. Here, we give a simple construction that deterministically eliminates 
all such anticrossings in polynomial time, space, and energy, for any Ising models with polynomial 
final gap. Thus, in order for adiabatic quantum optimization to require exponential time to solve 
any NP-complete problem, some quality other than this type of anticrossing must be unavoidable 
and necessitate exponentially long runtimes. 



INTRODUCTION 

The usefulness of adiabatic quantum optimization 
(AQO) for solving NP-complete problems, as originally 
proposed in [T], has been the subject of much debate in 
recent years. In AQO, the Hamiltonian of the quantum 
system begins as Hb, whose ground state is easy to pre- 
pare as the initial state of the system, and is evolved 
to a final Hamiltonian Hp, whose ground state is the 
optimal solution to an optimization problem. The run- 
ning time, t f , required to obtain a large amplitude of the 
ground state at the end of the evolution is proportional to 
fmin' where ^min is the minimum energy gap between the 
ground and first excited states, excluding excited states 
simply heading toward final degeneracy with the ground 
state, during the evolution. The complexity of AQO is 
therefore determined by the scaling of gmin with the prob- 
lem size. 

An example of an NP-hard [H] optimization problem 
well suited to being represented by Hp is that of find- 
ing the ground state of an Ising model with arbitrary 
couplings, Jij, and arbitrary local fields, hi: 

//p = 5]M^)+5]j.,aWap). (1) 

i i<j 

It is easily seen that this Hp is identical to the classical 
Ising model Hamiltonian, so its ground state is naturally 
identical to the ground state being sought (referred to 
here as the "global minimum" to avoid confusion where 
applicable). Hb is often chosen to be of the form 

HB = -^A,fT« A, >0, (2) 

i 

whose ground state is the in-phase, uniform superposi- 
tion of all Z-basis states. If representations of Hp and 
Hb can be found for a problem such that Qmin decreases 
only polynomially with the number of qubits, n, that 
problem can be solved in polynomial time using AQO. 
This, of course, assumes that the quantum system used 
to implement AQO is sufficiently isolated from a thermal 
environment, but effects of such an environment are not 
examined here. 
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FIG. 1: a) The two lowest energy eigenstates in an AQO spec- 
trum, which to low-order perturbation appear to cross, but 
actually anticross, causing a small Qmin, and b) three lowest 
eigenstates after changing the final ground state to be twofold 
degenerate. If the final ground state is degenerate, the cor- 
responding eigenstates (solid and dotted blue) will repel each 
other away from the end. If this degeneracy can be intro- 
duced without significantly affecting the excited state, the 
ground state energy (solid blue) will consequently be pushed 
away from that of the excited state (solid red) . 



The complexity of AQO for solving NP-complete prob- 
lems is a contentious issue, and has recently sparked in- 
tense debate. It was originally conjectured that the gap 
size may scale polynomially with problem size [T], but 
counterexamples were found [31 [3], which were subse- 
quently defeated by modifying the Hamiltonian [H [S]. 
Using perturbation expansion, Amin and Choi [B] showed 
that an eigenstate corresponding with a local minimum 
of Hp (anti)crossing with that of the global minimum 
near the end of the evolution, sometimes called a first 
order quantum phase transition, can result in an expo- 
nentially small 5min- In particular, when a local mini- 
mum has more low-energy, possibly degenerate, states in 
its vicinity than the global minimum, the eigenstate cor- 
responding to the local minimum may cross that of the 
global minimum. Using the same perturbation argument, 
Altshuler et al. [7] showed that for one representation 
of random exact cover instances, the probability of hav- 
ing problematic crossings increases with the system size, 
and claimed that these crossings are unavoidable. Oth- 
ers made similar observations using different techniques 
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Others argued that the analysis in [7j had not rigor- 
ously proven that the crossings are unavoidable, since 
problem structure that can impact the presence of such 
crossings was neglected [TUHT3]. Dickson and Amin [T3] 
showed that there is always some selection of Hb and Hp 
that guarantees no crossings of local and global minima, 
but it was left open whether or not an efficient method 
of selecting the Hamiltonian exists. 

Here, we present the first general method for eliminat- 
ing perturbative crossings, based on the effects of degen- 
eracy of the eigenstates of Hp, illustrated in Fig. [l] 

THE CONCEPT 

The main concept of the construction we present here 
is to force the ground state to diverge from all other states 
away from the end of the evolution by strategically intro- 
ducing degeneracy, as depicted in Fig. [l] and Fig. [5] ehm- 
inating any crossings that may have existed before apply- 
ing the construction. This can be done by adding extra 
qubits in such a way that the ground state becomes the 
most degenerate, the first excited state less degenerate, 
continuing this up to the highest excited state, which will 
be the least degenerate (possibly nondegenerate) . This 
degeneracy makes the eigenstates corresponding with the 
global minimum repel each other from the end more than 
all other eigenstates, so the ground state diverges from 
all other eigenstates as desired. This must be done with- 
out any knowledge of the global minimum, apart from 
that it is the lowest energy state of Hp. 

For simplicity, we will begin with the case of minimiz- 
ing arbitrarily-connected Ising models where all h,JE 
{—1,0,-1-1}, which is still NP-hard. This can be eas- 
ily generalized to arbitrary h ov J values (of polynomial 
precision) without resorting to expensive reduction map- 
pings onto {—1,0,-1-1}, as described later. Even though 
the only structure in a general, dense Ising model is that 
at most 2 spins participate in each term of its Hamilto- 
nian, this small amount of structure is enough to perform 
the required transformation. This does not apply to ran- 
dom energy models, studied in [TH [TS], since they have 
no applicable structure. 

The construction introduces extra qubits for each 
nonzero h ov J term of Hp such that when the h or 
J term is satisfied (i.e. at its negative value), the extra 
qubits corresponding to that term are at degeneracy (i.e. 
flipping them does not change the energy of the system), 
and when the term is unsatisfied, the extra qubits are not 
at degeneracy. This is shown for an h term in Fig. [2] and 
can be applied equivalently to a J term. Naturally, states 
with more terms satisfied will then be more degenerate 
than states with fewer terms satisfied, since states with 
more satisfied terms have more qubits at degeneracy. In 
the case where ft., J G {—1,0,-1-1}, the global minimum 



has the most terms satisfied, and will thus have the most 
qubits at degeneracy. As we will see, the slope (l^^'-order 
perturbation energy correction) of a state will be pro- 
portional to the number of qubits at degeneracy in that 
state, and so the ground state will have the largest neg- 
ative slope, diverging from all other states. 

THE CONSTRUCTION 

As mentioned above, we will begin by examining the 
case of Ising models with ft, J e {—1, 0, +1}. The Hamil- 
tonian of the original Ising model, M , is given as 

iyp = ^ft,aW+ J'^^^f^- (3) 

To introduce and amplify degeneracy of the ground 
state as desired, for each of the m nonzero terms in AI , 
we need to be able to add an extra qubit such that when 
the term is satisfied, the energy is unchanged regardless 
of whether the extra qubit is —1 or -t-1. However, when 
the term is unsatisfied, the energy must be unchanged 
when the extra qubit is —1, but higher by energy 26 > 
when the extra qubit is +1. This can be summarized as 





Added energy cost 


Term 


extra= 1 


extra=+l 


satisfied 








unsatisfied 





2b 



Note that if the energy didn't remain unchanged for 
some value of the extra qubit for both the satisfied and 
unsatisfied case, we would effectively be changing the en- 
ergy of the original states, which is not what we intend, so 
all three zeros are important. The added cost 2b serves 
only to break the degeneracy when the corresponding 
term is unsatisfied. 

To construct the new terms to appear in Hp, one may 
first observe that the above table is simply a multiple 
of the truth table for a boolean AND operator, which 
is equivalent to binary variable multiplication, so we 
need only to convert our two conditions ("unsatisfied" 
and "extra=-|-l") to binary variables and multiply by 2b. 
Since nonzero ft and J are limited to ±1, this gives us 
(term -I- l)/2 and (extra -I- l)/2, so the term added to 
the Hamiltonian for the extra qubit should be 

6(term+ 1) (extra + l)/2 (4) 

This is illustrated in Fig. [2] Now that we have a mech- 
anism for adding these extra qubits, for each hi ^ 0, add 
a > 1 extra qubits {ii,...,ia} € M^, with the following 
terms to the Hamiltonian. These are in exactly the same 
form as Q. 
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FIG. 2: Adding extra qubits in such a way that the ground state becomes degenerate, but the excited state remains nondegen- 
erate, causes the final ground states to repel each other. The ground state energy (sohd blue) then diverges from the excited 
state energy (solid red). Here, this is depicted for one original qubit (blue), and a) zero, b) one, and c) two extra qubits (green). 
A is the perturbation parameter used in the analysis below. For brevity, |t) <8> (|t) + li))/v^ is written as |t$), and similar for 
|t$$)- Note that the same method can be applied to a term +0"!"' instead of — trl"' simply by substituting one for the other 
in Hp. Likewise, the same method can be applied to a pair of original qubits coupled with a J term — tji^'o-i^' or -\-a'i^^ g'^zK 
again by substituting over — ai"'. 



^fe(;i,o-w + i)((T^) + i)/2 

k=l 

(5) 



Similarly, for each Jij ^ 0, add a extra qubits 
{iji, ija\ € ^i.i with the following terms to the Hamil- 
tonian. 



^ fo(^»j4'^^i'Vl)(a(^^'=^l)/2 (6) 

k=l 

It is important to note that while this introduces 3- 
local (ZZZ) terms, equivalent, albeit much less intuitive, 
2-local terms with the same property can be introduced 
by adding another extra qubit for each term. This is 
described in the appendix below. 

Since there are m nonzero terms in the original Ising 
model, M , we have added a total of am extra qubits. 
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each with an associated term coupUng it to M. 

Special consideration must be made when M may al- 
ready have single bit flip degeneracies. The above con- 
struction still works for sufficiently (polynomially) large 
a and b, but the analysis is made more complicated, so 
for simplicity of analysis in these cases, all qubits that 
may have a 1 bit flip degeneracy can be replaced by a 
strongly ferromagnetically coupled pair of qubits. Since 
the pair is stongly coupled to align with each other, they 
will act logically as the original single qubit, any 1 bit flip 
degeneracy with a 2 bit flip degeneracy (i.e. both qubits 
in the pair must flip to remain at the same energy). We 
are currently limiting the focus to J values of ±1, but 
this strong ferromagnetism can be expressed with multi- 
ple (identical) J = — 1 terms between the two qubits in 
the pair. Thus, the analysis below safely assumes that 
there are no 1 bit flip degeneracies in the original Ising 
model. 

Unfortunately, performing a full numerical simulation 
of this construction for relevant cases is infeasible, pri- 
marily due to the large number of added qubits neces- 
sary to apply it to Ising models sufficiently complex to 
have perturbative crossings. Thus, it is examined only 
analytically below. However, it is numerically feasible to 
demonstrate that changing the degeneracy, by adding ex- 
tra qubits in a similar manner, can be used to create and 
then eliminate these crossings. An example is illustrated 
in Fig. [3} 

ANALYSIS 

We will be performing a perturbation analysis on the 
above construction by expressing the time-dependent 
Hamiltonian, H, as 

H = Hp -f- \Hq 

i 

where A is initially oo, (though proportional to abm suf- 
fices), and is monotonically decreased toward zero over 
time. For now, we select all = 1. 

As mentioned above, the energies of the original states 
of AI (with the additional qubits in their —1 states) are 
unchanged by the addition of the new qubits, and all 
added states with higher energy are 1 bit fiip away from 
states with lower energy, and thus are not local minima. 

Supposing first that M has a nondegenerate state a, 
after applying the construction, a is now degenerate (un- 
less of course all terms were unsatisfied, which means that 
a is the highest energy state). In particular, upon intro- 
ducing small A > 0, the lowest eigenstate of the degener- 
acy, |a), is always a uniform superposition of a, since this 
uniquely minimizes {a \ Hb \q:), the coefficient of the 1^*- 
order term below. We wish to examine the perturbation 



expansion of this eigenstate \a), 

^^|a)(A)=i?fS+<^A + i?gA2 + ... (8) 

Because nonzero h, J are limited to ±1, we have that 

£'l'"| = — (# of terms satisfied in a) + (# unsat. in a) 
= —2{^ of terms satisfied in a) + m (9) 

Then, because |a) is a uniform superposition over the 
states corresponding with a, the f'^-order perturbative 
energy correction of \a) is 

- {a\HB\a) 

= - E A. 

i dcgcn. in a 

= —a{^ of terms satisfied in a) 

-I «,'-») 

= °j - constant. (10) 

This means that a final state a of lower energy than a 
final state f3, will have a lowest eigenstate \a) with larger 
negative slope than Thus, l'^*-order perturbation 

predicts that \a) and |/3) will diverge and not cross. 

In the case where a is degenerate in the original Ising 
model, (but with no 1 bit flip degeneracy) , it can be seen 
that I a) is a combination of uniform superpositions over 
each state's new degeneracy. For example, if a originally 
contained states k and I, \a) contains some combination 
of a uniform superposition over the states correspond- 
ing with k and a uniform superposition over the states 
corresponding with I. Since there are assumed to be no 
f^'-order degeneracies in the original a, this means that 

eI^I is the same as in the non-degenerate case. 
I") ° 

However, in order to prove that this can eliminate per- 
turbative crossings, one must examine the impact of this 
construction on higher orders of perturbation. Consider 
a single state a* in the eigenstate la). If a* is a local min- 
imum (a state with which a crossing could occur with the 
global minimum) , all states 1 bit flip from a* have higher 
energy, and therefore have at least one h or J unsatis- 
fied that is satisfied in a* . As decribed above, for any h 
or J satisfied in a* , each of the corresponding a added 
qubits contributes half of the states in the new superpo- 
sition I a*). This means that fiipping any 1 bit, i, from 
a*, which had incurred energy cost B^^^i, now incurs at 
least cost (1 + 2kb)Ba,.i for the portion of |a*) that has 
k added qubits in their states. Since the choice of a* 
was arbitrary, this applies to all of |a). For simplicity 
of analysis, we choose a — b — n? ^ oo, though this is 
extremely excessive in practice. We then find that the 
2nd Q]-(jg]- energy correction is bounded by 



5 



Visualization Hp Energy Spectrum 




, , 1- -10 

2 10 




FIG. 3: a) Starting from a simple 3-qubit Ising model without a perturbative (anti)crossing, b) the known local minimum, 
IH-i), can be made more degenerate (8-fold in this case) until a harsh crossing is created between the eigenstates corresponding 
with it and the global minimum, |ttt)- c) Now that there is a crossing, it can be eliminated again by increasing the degeneracy 
of the known global minimum to match or preferably exceed that of the local minimum. Although this does not apply the 
full construction (exact diagonalization is infeasible for the 42 qubits needed to eliminate the crossing), which does not require 
knowledge of the local or global minima, this example demonstrates the principle upon which the construction is based. (Note: 
For easier viewing, biases that exactly cancel are not depicted in the visualization.) 
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4 1") 



k=0 



E 



added i 
nondcgcn. in a 



2b 



1 (# of terms unsatisfied in a) 

ab 2b 



„4 ^ 4 V I") ^ J 



1 



1") 

constant. 



- Er'; + m 
4 V 1") 



(11) 



One may wonder why a = b was chosen instead of 
b ^ a, since the latter would make the 2"^^ order term 
approach zero, instead of favouring higher-energy final 
states as it currently does. That, however, could also 
make higher-order terms approach zero in a similar man- 
ner, dramatically increasing the radius of convergence 
of the series. By setting a = b and expanding a com- 
ponent of the perturbation, one finds that for q > 2, 
the qth-order term will have a component of magnitude 
(£:[°|/4)«-i (of alternat ing sign) and no components pro- 
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portional to a. For example, it is easily checked that 
the 3'''^ order term is dominated by (l/16)(i?|^°j)^. This 
means that the radius of convergence scales proportional 
to getting smaller as the problem size increases. 

Since the magnitude of the l'^' order term is propor- 
tional to a, and the value of a was chosen to be 8(n^), 
within the radius of convergence, 1*** order should domi- 
nate 0(n^£'|*-°j) orders of perturbation. 

Choosing such a large a also guarantees elimination of 
crossings if M had I''' order degeneracy of its local min- 
ima, since for any M with at least one non-zero term, 
^a^orig ^ After adding the extra qubits, the differ- 
ence between the 1^* order terms will still be il{n?), in 
favour of the global minimum. 

One may argue that adding qubits with such a large 
factor b could simply add a new crossing at a point where 
X ^ b and the original diagonal part of the Hamiltonian 
can be neglected, effectively giving a diagonal part of just 

a 

^^5(/i.aW + l)(aF + l)/2 + 

a 

E E^(-^'^-^i^^^i'^+i)('^i^""^+i)/2- (12) 

{i,j}eM k=l 

However, the ground state of this Hamiltonian is expo- 
nentially degenerate and one of these states is always 
trivially found. This is because there are no local min- 
ima: any nonzero term can be made zero by flipping the 
corresponding added qubit [i^. or ij^.) from +1 to —1. 

EXTENDING TO ARBITRARY COEFFICIENTS 

The above construction for Ising model terms with co- 
efficients in {—1, 0, -1-1} is most easily extended to integer 
coefficients. For a term with integer coefficient +k, sim- 
ply rewrite it as k identical terms with coefficient -1-1 and 
apply the construction as before on all k terms. The same 
splitting works with coefficients of — fc becoming k terms 
with coefficient —1. 

To extend this beyond integers to polynomial precision 
real numbers, one could simply rescale M to be polyno- 
mially approximated by integers, but a more reasonable 
construction is available from that 

- - E (13) 

i dcgon. in a 

One can simply round a coefficient down to the nearest 
integer below, apply the construction as before, then add 
a extra qubits as above whose Aj values are selected to 
be the remainder that was rounded off, instead of 1. This 
has negligible impact on higher order terms in that the 
radius of convergence is still proportional to 1 /-E'l^j j (un- 
less the coefficients were chosen to be scaled such that 



the bulk of weight lies in the rounded off components, 
which is easily fixed by scaling up the coefficients into a 
more reasonable range). 



CONCLUSION 

Above, we have presented a simple method of pos- 
ing the NP-hard problem of finding the ground state of 
an arbitrarily-connected Ising model with local fields as 
an adiabatic quantum optimization with no perturbative 
crossings between local and global minima. Thus, all NP- 
complete problems can be solved using adiabatic quan- 
tum optimization without encountering these crossings. 

It is critical to note that this does not prove that adia- 
batic quantum optimization can solve NP-complete prob- 
lems in polynomial time. However, it does mean that 
proving otherwise requires identifying some effect other 
than perturbative crossings that unavoidably results in 
exponentially long adiabatic runtimes. 
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APPENDIX: REDUCTION TO 2-LOCAL 

As mentioned above, the prescription above introduces 
a 3-local term when adding an extra qubit correspond- 
ing to J terms in the original Ising model. However, a 
2-local function with the necessary properties can be con- 
structed by adding another extra qubit. The properties 
to maintain are: 

1. all values of the function are > 

2. all configurations of qubits i and j that satisfy the 
associated Jij ^ term must have 2 configurations 
of the additional qubits, 1 bit flip apart, where the 
function is 

3. all configurations of qubits i and j that do not sat- 
isfy the associated Jij ^ term must have 1 config- 
uration of the additional qubits where the function 
is 

4. all other configurations must give values > 2b. 
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The following 4-qubit, 2-local function satisfies these 
properties for a positive J. {a^f' + l)/2 is abbreviated as 

X'i ^ etc . 

+6 (xj + + x^J^) (1 - 2xjj,,) + 8a;,j,* - 1 

+ (14) 

The first 2 lines inside the outer parentheses encode 
a +i7i*Vi'' Vi*""°^ term. Note that alternate encodings of 
this 3- local term where exactly 2 of i^, Xj, and Xij^ are 
flipped (e.g. replaced with 1 — etc.) do not maintain 
the needed properties of the overall function. The cost 
table of /(+) is as follows: 







/(+ 


/2b 








^ z — 


-1 


a^'"*^ = +1 


-1 


-1 


-1 







2 


-1 


-1 


+ 1 


2 




1 


-1 


+1 


-1 


1 







-1 


+1 


+ 1 


4 







+1 


-1 


-1 


1 







+1 


-1 


+ 1 


4 







+1 


+1 


-1 


4 







+1 


+1 


+ 1 


8 




1 



It's clearly visible that for some value of new qubit 
ijk*^ qubit ijk is at degeneracy when the original J term 
is satisfied, but when the J term is unsatisfied, neither 
qubit is at degeneracy. 

The following 4-qubit, 2-local function satisfies these 
properties for a negative J: 



^4 ((1 Xi)xj + XjXiji_ + Xijf,[l Xi)) 



b 
2 

+6{{l-Xi)+Xj + x^j^){l 



- 1 
(15) 



The first 2 lines inside the outer parentheses encode a 
term. Instead of flipping Xi (i.e. replac- 



ing with 1 — Xi), it also works to instead flip either Xj 
or Xij^, so long as it is only 1 of the 3. Note that an 
alternate encoding of this 3-local term where all 3 are 
flipped does not maintain the needed properties of the 
overall function. The cost table of /(_) as defined above 
is as follows: 





/(-) 


/2b 








u z — 


-1 


a^''"*^ = +1 


-1 


-1 


-1 


1 







-1 


-1 


+1 


4 







-1 


+1 


-1 


4 







-1 


+1 


+1 


8 




1 


+1 


-1 


-1 







2 


+1 


-1 


+1 


2 




1 


+1 


+1 


-1 


1 







+1 


+1 


+1 


4 








It's again visible that qubit ijk is at degeneracy when 
and only when the original J term is satisfied. 
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